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FRACTURE  TOUGHNESS  AND  CRACK  GROWTH  PROPERTIES  OF  ZEROPUR 


ABSTRACT 

The  fracture  toughness  and  crack  growth  parameters  of  Zerodurl,  a  low  expansion  glass  ceramic 
material,  were  determined.  The  fracture  toughness  was  determined  using  indentation  techniques 
and  was  found  to  be  0.9  MPa  m^^.  The  crack  growth  parameters  were  determined  using 
indented  biaxial  specimens  subjected  to  static  and  dynamic  loading  in  an  aqueous  environment. 
The  crack  growth  parameters  n  and  ln(B)  were  found  to  be  30.7  and  -6.837,  respectively.  The 
crack  growth  parameters  were  also  determined  using  indented  biaxial  specimens  subjected  to 
dynamic  loading  in  an  ambient  50%  relative  humidity  environment.  The  crack  growth 
parameters  n  and  ln(B)  at  50%  relative  humidity  were  found  to  be  59.3  and  -17.51,  respectively. 


INTRODUCTION 


Zerodurl  is  an  engineered  lithium  aluminosilicate  glass-ceramic  material  which  is  ideal  for  use 
as  telescope  mirrors  due  to  its  homogeneity  and  low  coefficient  of  thermal  expansion  [1].  The 
use  of  Zerodur  is  anticipated  in  the  Advanced  X-ray  Astrophysics  Facility  (AXAF)  and  the  Far 
Ultraviolet  Spectroscopic  Explorer  (FUSE).  Zierodur  is  also,  at  present,  the  pre-phase  B  baseline 
mirror  material  for  the  Orbiting  Solar  Laboratory  (OSL).  In  order  to  conform  to  fracture  control 
guidelines  [2]  the  components  made  from  the  Zerodur  material  must  be  either  proof  tested  and/or 
analytically  shown  to  be  capable  of  withstanding  the  anticipated  loading.  The  work  reported 
here  was  initiated  to  supply  the  necessary  material  dependent  crack  growth  data  to  enable 
analytical  predictions  of  component  lifetimes  to  satisfy  the  safe-life  design  criteria. 

To  demonstrate  safe-life,  the  initial  flaw  size,  the  flaw  growth  rate  and  the  critical  flaw  size  must 
be  known.  The  initial  flaw  size  is  dependent  on  the  manufacturing  processes  used  to  produce  the 
component.  The  strength  controlling  flaw  could  be  either  a  surface  flaw  created  during 
machining  the  component  or  volume  flaws  created  during  manufacturing  of  the  bulk  material. 

The  strength  controlling  flaws  in  glass  and  ceramic  components  are  generally  one  to  two  orders 
of  magnitude  smaller  in  size  than  those  that  can  be  located  with  conventional  non-destructive 
inspection  techniques.  Because  the  inherent  flaw  sizes  cannot  be  measured  direedy,  the  initial 
condition  of  the  component  cannot  be  directly  determined  (except  for  very  large  flaws).  The 
flaw  population  must  be  determined  indirectly  using  test  coupons  prepared  with  the  same  surface 
finish(es)  as  the  flight  component.  A  large  population  of  coupons  must  be  tested  to  accurately 
determine  the  spectrum  of  flaw  sizes  created  by  a  given  finislring  technique.  The  flaw  size  then 
used  as  the  initial  flaw  size  of  the  component  is  statistically  made  based  on  test  results  and  some 
acceptable  level  of  probability.  The  initial  flaw  size  must  be  determined  on  a  case  by  case  basis 
and  no  inference  as  to  initial  flaw  size  should  be  drawn  from  this  report. 

The  crack  growth  observed  in  glass  and  ceramic  materials  is  known  to  be  assisted  by  the 
presence  of  ambient  moisture.  To  calculate  the  amount  of  growth  a  flaw  will  under  go,  the  rate 
at  which  the  flaw  will  grow  in  the  service  environment  and  the  applied  stress  levels  must  be 
known. 
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The  fracture  toughness  of  the  material  must  also  be  known  to  demonstrate  safe-life.  The  frac^re 
toughness  of  the  material  will  dictate  how  much  stress  a  component  can  withstand  for  a  given 
flaw  size,  or  will  predict  the  flaw  size  that  will  induce  failure  for  a  given  stress. 

The  primary  purpose  of  the  testing  performed  for  this  report  was  to  determine  the  crack  growth 
rate  in  Zerodur.  The  fracture  toughness  has  also  been  determined  using  inderitation  techniques. 
While  specimens  were  tested  in  the  as-received  condition  to  quantify  the  initial  flaw  size 
distribution,  no  attempt  should  be  made  to  use  the  strength  of  the  as-received  specimens  as  a 
design  parameter. 


THEORETICAL  BACKGROUND 


Biaxial  Loading 

The  use  of  biaxial  loading  to  interrogate  glass  and  ceramic  materials  has  well  documented 
in  recent  years  [3-5].  The  primary  advantage  of  this  technique  is  the  specimen  cost.  A  cylmder 
with  the  desired  diameter  can  be  cote  drilled  from  bulk  material  and  the  disk  sj^cimens  sliced 
off  the  cylinder.  The  specimens  can  then  be  polished  to  the  desired  surface  finish.  Very  little 
attention  needs  to  be  paid  to  the  specimen  edges  because  the  creation  of  the  maximum  stress  at 
the  center  of  the  specimen  prevents  any  failures  from  edge  flaws,  as  is  observed  in  three  or  four 
point  bend  specimens. 

The  stress  on  the  tensile  surface  of  the  test  specimens  (a)  is  given  as  [3]: 


a  =  3P[2(H-v)ln(a/b)+(l-v)(a2-b2)/R2)]/4(K)t2  (1) 

where  P  is  the  applied  load  in  newtons,  v  is  the  Poisson’s  ratio,  t  is  the  specimen  thickness,  a  is 
the  support  ring  radius,  b  is  the  loading  ring  radius  and  R  is  the  specimen  radius.  The  specimen 
dimensions  are  in  millimeters. 

The  strength  data  at  each  test  condition  were  fitted  to  a  two  parameter  Weibull  relationship  [6]: 

In  In  [1/(1-F)]  =  m  ln[afi^ao]  (2) 


where  the  value  Of  is  used  to  denote  the  specimen  strength  at  failure,  m  and  Oo  are  the  Weibull 
modulus  and  Weibull  scale  parameter,  respectively.  The  Weibull  parameters  are  determined  by 
a  least  squares  linear  regression  of  Eq.  (2)  where  F  is  the  cumulative  failure  probability  defined 
as,  F=(i+0.5)/J  (i  is  the  specimen  rank  and  J  is  the  total  number  of  specimens  tested).  The 
Weibull  median  strength  is  the  fracture  strength  corresponding  to  a  failure  probability  of  0.5  and 
the  Weibull  modulus  represents  the  scatter  in  the  specimen  strengths. 

The  strength  of  large  components  can  be  predicted  based  on  small  specimen  strengths  by 
accounting  for  the  difference  in  size  or  area  stressed.  The  relationship  is  given  as; 

ai/02=(A2/Ai)l/«n  (3) 
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where  A  is  the  area  being  interrogated  and  m  is  the  WeibuU  modulus,  the  subscripts  1  and  2 
represent  the  large  component  and  test  specimens,  respectively.  Equation  (3)  accounts  only  for 
surface  flaws.  If  volume  flaws  are  the  strength  controlling  defects  then  the  volume  and  the  stress 
distribution  in  that  volume  must  be  taken  into  account. 

Indentation  Theory 

The  use  of  indenters  (Knoop  and  Vickers)  to  create  well  defined  flaw  systems  (Figure  1)  in  glass 
and  ceramic  components  has  been  well  documented  [7-14].  The  variation  in  strength  of 
specimens  with  indentation  flaws  is  much  lower  than  that  of  specimens  with  as-machined 
surfaces  and  is  thus  ideally  suited  for  die  introduction  of  flaws  for  crack  growth  testing.  The 
flaws  created  using  the  indentation  technique  are  on  the  scale  of  the  as  machined  flaws. 

The  response  of  brittle  materials  when  indented  with  a  Vickers  indenter  is  the  creation  of  a  well 
defined  crack  system  driven  by  contact  residual  stress.  The  contact  residual  stress  is  created  by  a 
strain  mismatch  between  the  plastically  deformed  material  adjacent  to  the  indentation  site  and 
the  elastically  deformed  surrounding  matrix.  During  the  indentation  event  two  crack  systems  are 
fomied;  radial  and  lateral.  The  radial  cracks  grow  peipendicular  to  the  surface  and  are  the  flaws 
from  which  fracture  occurs.  The  lateral  flaws  form  at  the  base  of  the  plastic  zone  and  grow 
parallel  to  the  surface. 


Figure  1.  Idealized  model  of  damage  created  by  indentation  with  Vickers  indenter.  The 
radial  crack  length  is  Co  and  the  lateral  crack  length  is  Cl.  The  daAened  area  beneath  the 
indentation  represents  the  plastic  zone  and  Oa-is  the  far  field  applied  stress  (after 
reference  [9]). 
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Fracture  Toughness 
Direct  Crack  Method 


The  stress  intensity  at  the  tip  of  the  radial  crack  created  during  indentation  is  determined  by 
modelling  the  indentation  stress  field  and  the  radial  crack  system  as  a  center  loaded  crack.  It  can 
be  shown  that  the  radial  crack  dimension  after  indentation  (cq)  is  [10]: 


Co  =  [XrP/Kic]2/3  (4) 


where  P  is  the  indenter  load,  Kic  is  the  critical  stress  intensity  factor  (fracture  toughness)  and  %x 
is  the  residual  stress  constant  that  is  dependent  on  indenter  geometry  and  material  properties. 
The  Xr  factor  has  been  calibrated  using  a  wide  range  of  ceramic  materials  and  was  found  to  be 
equal  to  0.016  (E/H)l/2  where  E  is  the  elastic  modulus  and  H  is  the  hardness.  Rearranging 
equation  (4)  yields  the  critical  stress  intensity  as  [10]: 


Kic  =  0.016[E/H]l/2[P/co3/2] 


(5) 


Note  that  the  [P/cq^/^]  term  is  a  material  constant.  This  value  should  be  constant  in  the  absence 
of  a  far  field  residual  stress  in  the  material  being  tested. 

Strength  Method 

The  stress  intensity  at  the  tip  of  a  radial  crack  may  be  found  by  summing  the  stress  intensity  due 
to  residual  and  applied  stresses.  The  post  indentation  strength  is  calculated  at  the  point  where 
crack  instability  occurs:  Ki=Kic  and  dKj/doO.  At  this  point  the  critical  stress  intensity  is  found 
to  be  [11]: 


Kic  =  Tl  [E/H]l/8  [afPl/3]3/4 


(6) 


where  tj  is  a  calibration  constant  that  also  inco)^orates  the  flaw  geometry,  assumed  to  be  an  ideal 
half  permy  shape.  The  T|  constant  has  been  calibrated  to  be  0.59  (±0. 12)  [1 1]. 

Crack  Growth  Parameters 

Static  Fatigue 

Crack  growth  of  glass  and  ceramic  materials  occurs  when  statically  loaded  in  an  aggressive 
environment.  The  crack  growth  can  be  expressed  as  a  power  function  of  the  stress  intensity 
factor  (Ki)  and  is  given  as  [14-17]: 


V  =  VO  (Ki/Kic)i^ 

where  v  is  the  crack  velocity,  and  vq  &  n  are  material  and  environmental  crack  growth 
constants.  From  the  power  law  relationship  an  expression  for  the  time  to  failure  (tf)  can  be 
derived  in  terms  of  the  applied  stress  (<Ja)>  is  given  as: 
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tf  =  B  ai"‘2  aa‘” 

B  =  2/[voY2  (n-2)Kic(-2)] 


(8) 

(8a) 


where  ai  is  the  strength  in  the  absence  of  subcritical  crack  growth  or  inert  strength  and  Y  is  a 
flaw  geometry  factor  taken  to  be  rt  for  surface  flaws. 

The  crack  growth  parameters  n  and  B  are  found  using  the  median  time  to  failure  and 
corresponding  applied  stress.  The  use  of  median  values  are  convenient  as  they  represent  the  time 
to  failure  at  a  50%  failure  probability.  A  least  squares  linear  regression  analysis  is  performed  on 
the  natural  logs  of  the  median  time  to  failure  and  the  corresponding  applied  stress.  The  slope  is 
the  n  parameter  and  the  intercept  is  [B  Oi“‘^]. 

Dynamic  Fatigue 

D)mamic  fatigue  of  glass  and  ceramic  materials  refers  to  the  decrease  in  the  observed  strength  as 
a  function  of  decreasing  stressing  rate.  The  relationship  between  failure  stress  (Of)  and  stressing 
rate  (o)  is  given  as: 


Q^+l  =  B  (n+1)  oi*'"^  O  (9) 

The  crack  growth  parameters  B  and  n  are  determined  as  described  above  for  static  fatigue,  using 
the  median  failure  stress  and  the  associated  stressing  rate. 


Fatigue  of  Indentation  Flaws 

.It  has  been  shown  that  the  fatigue  parameters  determined  using  indentation  flaws  differ  from 
those  determined  when  the  residual  stress  field  associated  with  the  indentation  event  is  removed 
by  annealing.'  The  relationship  between  the  apparent  fatigue  exponent  (n')  and  actual  fatigue 
exponent  (n)  is  given  as  [17]: 


n  =  1.31  n' 


(10) 


The  primed  notation  designates  the  value  determined  in  the  presence  of  the  residual  stress.  The 
relationship  between  the  B  and  B '  values  is  not  as  straight  forward  as  that  of  n.  The  true  value 
of  the  fatigue  parameter  vq  as  a  function  of  indentation  load  is  given  as  [17]: 


Vo  =  [2.84  n '  0.462  <jj^n '  Cj^]  A'  (11) 

X'  =  B '  (n'  +1)  -^)  (1  la) 


where  am  is  the  inert  strength  for  a  given  indentation  load  and  c^^  is  the  corresponding  radial 
crack  size  at  failure  and  taken  to  be  2.5  times  larger  than  Co-  The  value  for  am  is  taken  as  die 
inert  strength  in  equations  (8)  and  (9).  The  trae  B  parameter  can  be  determined  by  substituting 
values  obtained  in  equations  (10)  and  (11)  into  equation  (8a). 
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TEST  MATERIAL 


Zerodur  is  an  opaque,  engineered  glass-ceramic  material  based  on  a  lithium-aluminosilicate 
formulation  with  the  crystalline  phase  consisting  of  metastable  solid  solutions  of  high  quartz 
stmcrnre.  Zerodur  is  70  to  78%  by  weight  crystalline  phase  with  crystals  generally  50  to  55 
nanometers  in  size.  The  material  composition  is  given  in  Table  1. 


Table  1.  Zerodur  Composition*  Table  2.  Zerodur  Mechanical  Properties* 


Composition 

Modulus  of  Elasticity  (E,GPa) 

91.0 

wt% 

mol% 

Poisson’s  Ratio  (v) 

0.24 

Density  (p,g/cm^) 

2.53 

Si02 

55.4 

63.9 

Hardness,  Knoop  (H,  GPa) 

6.18 

AI2O3 

25.4 

17.2 

Strength  (MPa) 

90.0** 

Li20 

3.7 

8.5 

Na20 

0.2 

0.2 

*  Manufacturers  data. 

K20 

0.6 

0.5 

**  0,5  cm2  test  area,  loaded  at  2  MPa/sec,  No.  600 

MgO 

1.0 

1.7 

loose  grain  finish,  5%  failure  probability. 

ZnO 

1.6 

1.3 

P2O5 

7.2 

3.5 

Ti02 

2.3 

2.0 

Zr02 

1.8 

1.0 

As03 

0.6 

0.2 

*  Ref.  [19] 


The  manufacture  of  the  Zerodur  material  requires  tight  control  of  a  nucleation  and  crystallization 
process  which  creates  the  material’s  microstructure.  Annealing  of  sodalime  glass  specimens 
with  indentation  flaws  has  been  shown  to  remove  the  localized  residual  stress  created  during  the 
indentation  event  [15].  Annealing  of  the  Zerodur  specimens  to  remove  the  indentation  residual 
stress  field  might  have  altered  this  microstructure,  so  was  not  performed. 

Zerodur  has  very  good  thermal  stability.  Various  thermal  expansion  classes  are  available  with 
the  most  stable  being  expansion  class  1  which  has  a  thermal  expansion  coefficient  of  0  ± 
0.05ppm/K  for  0  to  50oC. 

Zerodur  has  good  resistance  to  chemical  attack.  Zerodur  resists  attack  by  most  acids  (with 
exception  of  hydrofluoric  acid),  alkalis  and  salts  at  room  temperature.  Sulfuric  acid  does  attack 
Zerodur  at  higher  temperatures. 


TEST  SPECIMENS  AND  PROCEDURE 

Zerodur  disks  (38.1  mm  dia  x  1.78  mm)  were  used  in  this  study.  The  faces  were  polished  to  a 
commercial  grade  finish  by  the  manufacturer!.  All  specunens  were  randomized  prior  to  testing 
to  avoid  any  systematic  error  due  to  manufacturing  processes. 
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Indentation 


A  Vickers  diamond  indenter  was  used  to  introduce  indentations  into  the  Zerodur  specimens  at 
loads  ranging  from  300  to  2000  grams.  The  hardness  tester^  used  to  introduce  indentations  into 
the  specimens  had  a  constant  loading  rate  of  6  mm/min  and  remained  in  contact  for  20  seconds. 
Crack  lengths  were  measured  one  hour  after  the  indentation  event  on  a  Nikon  Epiphot 
microscope  at  1000  times  magnification. 

The  indents  were  placed  in  the  center  of  the  tensile  side  of  the  specimens.  The  center  of  the 
specimen  was  located  by  placing  the  specimen  on  a  template  and  placing  a  small  fiducial  "dot" 
in  the  center  of  the  compressive  face.  The  dot  was  a  dry  transfer  lettering  bullet.  The  center  of 
the  specimen  was  then  located  directly  under  the  indenter.  This  fiducial  dot  also  insured  that  the 
center  of  the  specimen,  and  indentation,  were  located  at  the  center  of  the  loading  ring.  The  test 
specimens  were  allowed  to  sit  in  ambient  conditions  for  one  hour  prior  to  strength  testing. 

Dynamic  Testing 

A  biaxial  test  fixture  mounted  on  a  servo  hydraulic  test  machine^  was  used  for  all  dynamic 
strength  testing.  All  specimens  were  tested  with  a  ring  on  ring  loading  configuration  (Figure  2). 
The  loading  and  support  rings  were  5.8  and  25.4  mm  diameter  thrust  bearings,  respectively. 

The  use  of  the  thrust  bearings  to  support  the  specimens  was  preferred  over  machined  rings  for 
several  reasons.  The  use  of  thmst  bearings  was  felt  to  lessen  the  amount  of  friction  on  the 
specimen  surface,  there  is  a  high  degree  of  dimensional  consistency  between  bearings  and  the 
cost  of  the  thmst  bearings  was  much  less  than  the  cost  of  machming  several  rings  The  test 
fixtures  used  were  easily  retrofitted  using  the  thmst  bearing  races.  The  smaller  thmst  bearings 
were  used  as  loading  rings  to  avoid  any  effect  that  might  arise  at  the  periphery  of  a  ground  ball 
and  to  provide  dimensional  consistency  between  loading  fixtures. 


Lxjading  Shaft 

1/4**  diameter  ball 

Bearing  Race 
Loading  ThniSl  Bearing 

Test  Specimen 

Support  Tliru St  Bearing 

Bearing  Race 

Cotton  Wick 
Moisture  Reservoir 
Plastic  Tubing 


Figure  2.  Schematic  of  test  configuration. 


2  Tukon  Tester  Model  FB,  Wilson  Instruments,  Binghamton,  NY 

3  InstTon  1350,  Instron  Corp.,  Canton,  MA 
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The  inert  strength  of  the  Zerodur  specimens  was  detennined  for  both  as-received  and  indented 
sf>ecimens.  The  indentation  loads  were  300,  500,  1000  and  2000  grams.  A  minimum  of  10 
specimens  was  tested  at  each  condition.  The  inert  strength  specimens  were  tested  in  air  at  a 
loading  rate  of  2  mUlimeters/second  which  resulted  in  stressing  rates  of  approximately  2000 
MPa/sec.  This  loading  rate  was  felt  to  be  sufficient  to  limit  the  environmentdly  assisted  crack 
growth. 

The  load  was  recorded  using  a  data  acquisition  and  plotting  system^.  The  load  was  monitored 
with  a  piezoelectric  load  cell5  at  loading  rates  greater  than  0.2  mm/sec  and  with  a  standard  strain 
gauged  load  cell  at  loading  rates  less  than  0.2  mm/sec.  The  slowest  loading  rate  used  was  5 
mm^. 

The  dynamic  fatigue  parameters  were  determined  for  a  water  environment  using  specimens  with 
500  and  1000  gram  Vickers  indentations.  The  water  environment  was  created  by  placing  a  drop 
of  deionized  water  over  the  indentation  prior  to  testing.  The  dynamic  fatigue  parameter  was  also 
determined  for  50%  relative  humidity  (RH)  using  specimens  with  1000  gram  indentations.  The 
50%  RH  was  provided  by  the  ambient  lab  environment. 

Static  Testing 

The  static  fatigue  testing  reported  here  used  only  two  fixtures.  The  fixtures  used  for  static 
fatigue  testing  were  identical  to  that  used  for  the  dynamic  fatigue  testing  but  had  a  weight  pan 
attached  to  the  loading  shaft.  It  was  confirmed  that  the  applied  stress  on  the  tensile  surface 
agreed  with  theory  using  a  strain  gauged  Zerodur  specimen.  The  agreement  with  theory  was 
within  1  percent  for  both  fixtures. 

A  target  stress  was  selected  and  the  desired  weight  calculated  using  Eq.  (1)  and  the  individual 
specimen  thickness.  The  load  was  applied  with  lead  blocks  and  lead  shot.  The  lead  shot  was 
weighed  out  to  the  precise  load  required  to  induce  the  target  stress.  The  load  was  applied  to  the 
specimen  via  a  weight  pan.  When  specimen  failure  occurred  the  weight  pan  tripped  a 
microswitch  and  the  failure  time  was  automatically  recorded  by  a  personal  computer^. 

The  static  fatigue  specimens  were  tested  in  a  deionized  water  environment.  The  water  was 
wicked  to  the  specimen  surface  from  a  reservoir  located  directly  beneath  the  tensile  surface.  The 
reservoir  was  replenished  as  required  via  plastic  tubing. 


RESULTS 


Hardness 


The  indentation  diagonal  sizes  are  presented  in  Table  3.  The  hardness,  calculated  as  twice  the 
applied  load  divided  by  the  indentation  diagonal  squared  (the  factor  of  two  represents  the 
Vickers  indenter  geometry),  was  found  to  decrease  with  increasing  indenter  load.  This  variation 
of  hardness  with  increasing  load  has  been  previously  noted  with  sodalime  glass  [21].  The  mean 
hardness  value  was  found  to  be  6.37  GPa.  This  value  agrees  with  the  hardness  value  provided  by 
the  manufacturer  (6.2  GPa). 


Hewlett  Packard  Measurement  Plotting  System  (HP-7090A),  Hewlett-Packard  Co.,  Rockville  MD 
Kistler  Instrument  Corp.,  Amherst  NY 
COMPAQ,  COMPAQ  Computer  Corp., 
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Table  3.  Indentation  Dimensional  Measurements 


Indentation 

Load 

(grams) 

Indentation 

Diagonal 

(um) 

Hardness 

(GPa) 

Crack 

Length 

(co,um) 

Fracture* 

Toughness 

(MPaml/2) 

300 

29.2(1.1) 

6.89 

32.1(1.8) 

0.941 

500 

39.2(0.9) 

6.39 

46.3(1.6) 

0.940 

1000 

56.3(1.1) 

6.20 

74.7(3.2) 

0.931 

2000 

80.8(2.9) 

6.01 

129.4(9.8) 

0.830 

*  Calculated  using  equation  (5) 

Numbers  in  parenthesis  are  standard  deviations. 


Inert  Strength 

The  inert  strength  data  is  presented  in  Table  4.  The  as-received  strength  value  determined  in  this 
study  can  be  compared  to  previously  reported  Zerodur  strength  (90  MPa,  Table  2)  by  first 
calculating  the  strength  at  a  5%  failure  probability  using  Eq.  (2)  and  then  scaling  the  5%  failure 
probability  strength  to  the  appropriate  test  area  using  Eq.  (3).  Taking  ao=197.9  MPa  and  m=3.2 
(Table  4)  the  as-received  strength  at  a  5%  failure  probability  is  78.2  MPa.  Then  scaling  from  the 
0.264  cm2  area  used  in  this  study  is  to  the  0.5cm2  area  previously  used  the  strength  is  64.1  MPa. 
The  lower  strength  determined  from  this  study  indicates  that  the  surface  finish  was  slightly  better 
on  the  previously  tested  specimens. 


Table  4.  Inert  Strength 


WeibuU 

Parameters 


Indent 

Load 

(grams) 

Specimens 

Tested 

Median 

Strength 

(MPa) 

Mean 

Strength 

(MPa) 

KiC* 

(MPa-m^^) 

CTo 

(MPa) 

m 

As-received 

18 

154.9 

175.2(69.) 

197.9 

3.2 

300 

10 

105.9 

107.2(5.6) 

1.101 

109.9 

21.5 

500 

10 

86.6 

87.7(2.5) 

1.092 

89.1 

36.2 

1000 

20 

72.0 

71.0(5.3) 

1.141 

73.5 

15.6 

2000 

10 

57.8 

57.9(2.7) 

1.156 

59.2 

24.3 

*  Calculated  using  equation  (6) 

The  numbers  in  parenthesis  represent  1  standard  deviation 


The  mean  and  standard  deviation  data  are  plotted  in  Figure  3  on  a  In-ln  plot.  The  line  through 
the  data  in  Figure  3  has  a  slope  of  -1/3  with  the  [afP^^]  term  taken  to  be  152.7  MPaNl/3. 
Indentation  theory  predicts  that  the  strength  of  an  indented  material,  in  the  absence  of  any  far 
field  residual  stress,  will  vary  as  the  indentation  load  to  the  1/3  power.  The  test  results  are  in 
good  agreement  with  theory. 
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Figure  3,  Inert  specimen  strength  of  the  Zerodur  as  a  function  of  the  indentation  load. 
The  shaded  area  to  the  left  represents  the  mean  and  standard  deviation  of  the  ^-received 
specimen  strength.  The  line  through  the  data  has  a  theoretical  slope  of  -1/3  with  the 
constant  (aP'^)  value  taken  to  be  152.7  MpaN*/3. 


Fracture  Toughness 

The  fracture  toughness  determined  using  the  direct  crack  method  and  the  strength  method  were 
calculated  with  equations  (5)  and  (6),  respectively,  and  are  presented  in  Tables  3  and  4.  The 
average  fracture  toughness  determined  by  the  direct  crack  method  was  0.911  MPa-m^^  and  the 
average  fracture  toughness  determined  from  the  strength  method  was  1.123  MPa  m^^.  These 
fracture  toughness  values  are  greater  than  previously  determined  using  a  large  crack 
measurement  technique  (0.844  MPa-m^^)  [22]. 

Fatigue  Testing 

The  static  fatigue  median  value  data  and  associated  standard  deviations  are  plotted  in  Figure  4. 
The  fatigue  parameters  calculated  from  the  median  value  analysis  of  the  static  fatigue  d^  arc 
listed  in  Table  5.  The  data  shown  in  Figure  4  at  the  25  MPa  r^rplied  stress  was  not  used  to 
calculate  the  fatigue  parameters.  Testing  at  this  stress  level  is,  at  the  time  of  this  writing,  still 
under  way  (three  specimens  failed  prior  to  predicted  median  time  to  failure  and  two,  currently 
being  tested,  have  exceeded  median  time  to  failure). 
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Figure  4.  Static  fati^e  data  for  Zerodur  with  9.8  newton  indentations  in  a  water 
environment.  The  line  through  the  data  was  best  fit  to  the  median  values  (triangles).  The 
error  bars  represent  one  standard  deviation  and  the  crosses  represent  individual  failures. 
The  specimens  at  25  Mpa  were  not  used  in  the  regression  analysis. 


The  results  of  dynamic  fatigue  testing  using  median  value  analysis  are  presented  in  Figures  5  and 
6  and  the  calculated  fatigue  parameters  are  listed  in  Table  5.  The  fatigue  parameter  value  listed 
for  the  1000  gram  indentations  tested  in  50%  relative  humidity  was  calculated  only  from 
specimen  strengths  when  the  rate  of  crack  growth  was  dominated  by  the  ambient  inoisture. 
These  data  are  at  the  four  lower  stressing  rates  shown  in  Figure  6.  The  crack  growth  at  the  three 
higher  stressing  rates  in  Figure  6  is  not  completely  controlled  by  ambient  moisture. 


Table  5.  Apparent  Fatigue  Parameters  Determined  from  Median  Values 


Indentation 

Load 

(grams) 

Loading 

Environment 

Fatigue  Parameters 
lii(B')  n' 

500 

Dynamic 

Water 

-2.192 

23.6 

1000 

Dynamic 

Water 

-3.631 

23.5 

1000 

Static 

Water 

-8.252 

30.6 

1000 

Dynamic 

50%  RH 

-13.885 

45.2 
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DISCUSSION 


While  the  as-received  strength  of  the  specimens  tested  in  this  study  should  not  be  used  as  the  as- 
received  strength  of  flight  components,  it  is  instructive  to  observe  the  strength  distribution  of  the 
specimens  tested  in  the  as-received  condition.  The  Zerodur  specimens  used  in  this  study  were 
polished  to  a  commercial  grade  finish.  The  Weibull  plots  for  the  inert  strengths  are  shown  in 
Figure  7.  The  distribution  of  the  as-received  strengths  appears  to  be  that  of  a  bimodal 
distribution.  Bimodal  distributions  generally  occur  when  the  strength  controlling  flaws  come 
from  two  distinct  flaw  populations.  These  populations  might  be  inherent  voids  or  bubbles 
created  during  the  pouring  of  the  Zerodur  and  line  or  point  flaws  created  by  machining  or 
handling  of  the  specimens.  The  distinct  flaw  populations  could  also  arise  due  to  specimen 
preparation  preformed  by  different  machinists,  on  different  machines  or  on  different  days.  The 
fracture  origin  could  only  be  located  in  the  weakest  of  specimens  and  that  defect  appears  to  be  a 
linear  flaw. 


Figure  7.  Weibull  distribution  of  inert  strength  data. 


t 


13 


The  variation  of  fracture  toughness  values  obtained  from  different  techmques  (Table  6) 
illustrates  the  uncertainty  in  any  value  obtained  using  a  single  technique.  As  demonstrated  by 
the  fractme  toughness  values  determined  for  sodalime  glass,  the  indentation  strength  technique 
gives  an  approximate  fracture  toughness  value  but  is  generally  found  to  overestimate  the 
generally  accepted  fracture  toughness  values  for  very  low  Kjc  glass  materials  (<lMPa  m^^). 
The  over  estimation  of  the  fracture  toughness  of  glass  materials  by  the  strength  method  may  be 
due  to  the  range  over  which  the  calibration  constant  is  obtained.  The  calibration  constant  for  the 
strength  method  is  determined  using  materials  ranging  from  0.7  MPa  m^/^  to  10  MPa  m*/^  [11]. 
The  standard  deviation  associated  with  the  calibration  constant  is  20%  of  the  mean  value.  Due 
the  consistent  overestimation  of  the  fracture  toughness  of  low  Kic  materials  when  using  the 
indentation  strength  technique,  it  would  be  pradent  to  use  a  twenty  percent  lower  calibration 
constant  to  calculate  the  fracture  toughness.  Using  a  calibration  constant  of  0.47  (in  place  of  the 
mean  value,  0.59)  in  Eq.(5)  the  ^rodur  fracture  toughness  determined  from  the  strength 
technique  is  0.90  MPa  m^/^,  which  is  in  good  agreement  with  that  determined  by  the  direct  crack 
method. 

The  fracture  toughness  calculated  using  the  direct  crack  method  was  slightly  greater  than  that 
using  a  large  flaw  technique.  The  direct  crack  method  usually  is  in  good  agreement  with  the 
large  flaw  technique.  The  fracture  toughness  of  sodalime  glass  determined  by  the  direct  crack 
method,  using  the  same  equipment  and  procedures  as  was  used  to  test  the  Zerodur,  was  found  to 
agree  with  values  determined  at  other  laboratories.  This  discrepancy  may  be  due  to  the 
specimens  originating  from  different  material  lots.  To  insure  the  component  design  is 
conservative  the  lower  toughness  value  (0.844  MPa  m^^)  should  be  used  for  design  purposes. 


Table  6.  Comparison  of  Fracture  Toughness  Values(MPa  m'^) 


Zerodur 

Soda-Lime  Glass 

Ref  19 

GSFC(23] 

Direct  Crack 

0.911 

0.70 

0.67 

Strength  Method 

1.123 

0.94 

0.94 

Large  Crack  Method 

0.844(19] 

0.75 

Upon  initial  examination  of  the  apparent  fatigue  parameters  in  Table  5,  there  appears  to  be  poor 
agreement  between  fatigue  parameters  determined  from  static  and  dynamic  testing.  The 
discrepancy  is  eliminated  when  the  adjusted  dynamic  fatigue  parameters  are  compared  to  the 
unadjusted  static  fatigue  parameters  (Table  7).  The  difference  observed  in  the  fati^e  parameters 
may  be  explained  by  the  relaxation  of  the  localized  residual  stress  field  associated  with  the 
indentation  event  due  to  the  lateral  crack  growth  during  the  early  stages  of  static  loading. 

Recent  studies  have  shown  that  the  strength  of  sodalime  glass  specimens  with  indentation  flaws 
increases  with  time  after  the  indentation  event  [24].  It  was  observed  that  both  the  radial  and 
lateral  cracks  grew  during  exposure  to  aqueous  environments,  but  the  strength  increase  has  been 
attributed  to  the  lateral  crack  growth  relieving  the  localized  residual  stress.  While  the  strength  of 
the  specimens  tested  in  reference  [24]  did  increase,  the  maximum  strengths  were  less  than  that  of 
the  indented  and  annealed  specimens.  The  strength  difference  between  the  aged  and  aimealed 
specimens  is  attributed  to  the  growth  of  the  radial  cracks  during  aging. 

If  it  is  assumed  that  the  sustained  static  loading  allowed  the  lateral  cracks  to  grow  and  that  the 
localized  residual  stress  was  relieved  then  the  fatigue  parameters  determined  from  the  static 
fatigue  testing  are  the  true  fatigue  parameters  and  do  not  need  to  be  adjusted.  The  adjusted 
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parameters  from  the  dynamic  data  agree  quite  well  with  the  non  adjusted  parameters  determined 
from  the  static  fatigue  testing  (Table  7).  The  slightly  lower  B  value  for  the  static  fatigue  data 
may  be  attributed  to  the  slightly  different  flaw  geometry/stress  field  resulting  from  Ae  growth  of 
the  lateral  flaws;  i.e.  the  determination  of  the  B  parameter  is  dependent  on  the  inert  strength 
value. 


Table  7.  Fatigue  Parameters  in  the  Absence  of  Residual  Stress 


Indentation 

Loading 

Environment 

Fatigue  Parameters 

Load 

(grams) 

ln(B)* 

500 

Dynamic 

Water 

-5.347 

30.9 

1000 

Dynamic 

Water 

-6.913 

30.8 

1000 

Static** 

Water 

-8.252 

30.6 

1000 

Dynamic 

50%  RH 

-17.51 

59.2 

*  True  B  parameter  calculated  using  equations  (8)  &  (1 1). 
Static  parameters  not  adjusted 


It  seems  justified  at  this  point  to  draw  the  conclusion  that  the  results  obtained  using  static  and 
dynamic  testing  do  agree.  The  average  fatigue  parameters  for  Zerodur  in  the  absence  of  residual 
stresses  are  n=30.7  and  the  ln(B)  =  -6.837.  For  design  purposes  the  B  value  determined  from  the 
static  fatigue  testing  could  be  used  to  insure  a  conservative  design.  To  further  insure  that  the 
design  is  conservative  the  fatigue  parameters  determined  in  the  presence  of  the  indentation 
localized  residual  stress  field  (Table  5)  could  be  used  for  life  time  predictions.  The  use  of  the 
latter  parameters  could  allow  for  the  presence  of  residual  stresses  created  during  machining  of 
the  component,  which  would  increase  the  crack  growth  rate. 

While  it  is  necessary  to  insure  the  component  wiU  survive  launch,  excessive  conservatism  will 
add  undesired  weight  to  the  component.  The  use  of  crack  growth  parameters  deterr^ed  in  a 
water  environment  to  predict  component  lifetimes  adds  a  greater  degree  of  conservatism.  It  is 
clear  from  the  data  presented  in  Tables  5  and  7  that  the  rate  of  crack  growth  in  ambient  moisture 
is  much  less  than  that  in  an  aqueous  environment.  The  values  determined  for  ln(B' )  and  n'  at 
50%  RH  were  -13.89  and  45.2,  respectively.  When  these  values  are  adjusted  to  reflect  the  true 
values,  ln(B)  and  n  values  are  -17.51  and  59.2,  respectively.  These  values  are  in  fair  agreement 
with  values  determined  in  reference  [22]  which  used  "...unilaterally  toothed,  tensile  samples 
loaded  on  one  axis..."  in  an  ambient  environment  of  40  to  60%  RH.  In  that  reference  the  ln(B) 
an  n  parameters  were  found  to  be  -11.49  and  51.7,  respectively.  The  difference  in  fatigue 
parameters  may  again  be  due  to  batch  to  batch  material  variations. 

It  is  clear  from  the  above  results  that  the  actual  component  environment  should  be  taken  in  to 
account  during  component  design.  Figure  8  shows  curves  generated  using  equation  (8),  an 
assumed  proof  stress  level  of  70  MPa  and  the  parameters  listed  in  Table  7.  The  initial  or  inert 
strength  in  equation  (8)  can  be  replaced  with  the  proof  stress  and  the  time  calculated  is  a 
minimum  time  to  frdliu-e.  The  70  MPa  stress  was  arbitrarily  selected  and  corresponds  to  a  failure 
probability  during  proof  stressing  of  3.7%.  It  can  be  seen  in  Figme  8  that  the  variation  of  the  B 
parameter  does  not  significantly  affect  the  lifetime  analysis  while  the  component  environment 
does  severely  affect  the  component  lifetime. 
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Figure  8.  Median  specimen  lifetime  calculated  using  fatigue  parameters  listed  in  Table  8 
and  in  reference  [22].  Inert  strength  was  assumed  to  be  70  Mpa. 


CONCLUSIONS 


The  fracture  toughness  determined  from  this  study  was  found  to  be  approximately  0.9  MPa  m*/^. 
A  previous  study  [19]  using  a  large  crack  propagation  technique  found  the  fracture  toughness  to 
be  0.844  MPa  m^^.  It  is  felt  the  lesser  value  should  be  used  to  insure  a  conservative  design. 

The  average  true  crack  growth  parameters  n  and  ln(B)  for  Zerodur  in  100%  relative  humidity 
were  determined  to  be  30.7  and  -6.837,  respectively.  The  tree  crack  growth  parameters  n  and 
ln(B)  for  Zerodur  in  50%  relative  humidity  were  determined  to  be  59.3  and  -17.51,  respectively. 
In  the  presence  of  a  residual  stress  these  v^ues  are  considerably  less. 
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ADDENDUM 


Homologous  Ratio  Analysis 

A  material’s  crack  growth  constants  are  generally  determined  using  the  median  values,  since 
these  values  are  at  equal  failure  probabilities  (0.5).  The  median  values  analysis,  while  easy  to 
apply,  can  result  in  large  uncertainties  in  the  B  and  n  values.  An  alternate  method  to  analyze  the 
data  is  based  on  a  homologous  stress  ratio  (<Jhs).  given  as  [20]: 


CtHS  =  <ya/Oi 


(11) 


The  homologous  stress  is  ranked  based  on  the  inert  strength  ranking  and  associated  with  the 
correspondingly  ranked  time  to  failure  for  a  given  applied  stress.  Equation  (8)  is  rewritten  as: 

InCtfCTi^)  =  InB  -  n  In  ojis  (12) 


The  dynamic  fatigue  can  be  similarly  analyzed  by  defining  the  homologous  stress  ratio  to  be: 


ctHD  =  cy£(<yi 


(13) 


and  equation  (9)  can  be  rewritten  as: 

In  <thd  =  (n+l)-l[ln  B  +  In(n+1)  +  ln(a/(Ji3)] 


(14) 


Static  Fatigue 

The  static  fatigue  data  analyzed  using  the  homologous  ratio  analysis  technique  are  presented  in 
Figure  9.  The  fatigue  parameters  determined  at  each  applied  stress  and  that  determined  using  all 
the  data  are  listed  in  Table  8.  When  the  ranked  data  had  unequal  populations,  data  from  the 
larger  set  were  randomly  selected  and  removed  from  the  population,  the  population  was  then  re¬ 
ranked. 

The  cause  for  the  variation  in  the  fatigue  parameters  obtained  using  the  homologous  ratio  to 
analyze  the  static  fatigue  data  is  not  certain.  The  homologous  static  fatigue  parameters 
determined  using  all  stress  levels  were  close  to  the  parameters  calculated  using  the  median 
values.  Two  individual  stress  levels  were  in  excellent  agreement  while  the  other  two  deviated 
somewhat.  This  technique  uses  all  the  data  but  it  is  not  known  if  there  is  a  minimum  number  of 
data  points  required  to  insure  convergence  to  the  correct  fatigue  parameters. 
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Table  8.  Fatigue  Parameters  Determined  from 
Homologous  Ratio  Analysis  of  Static  Fatigue  Data 


Applied 

Stress 

(MPa) 

Fatigue  Parameters 
lii(BO  n' 

28 

-10.730 

35.1 

30 

-8.736 

30.7 

32 

-14.190 

37.9 

36 

-7.645 

30.9 

Combined 

-8.579 

31.5 

The  homologous  ratio  analysis  of  the  dynamic  fatigue  data  at  individual  stress  rates  did  not  yield 
realistic  values  for  the  crack  growth  parameters.  TTie  values  for  the  n  parameter  calculated  from 
the  500  gram  indentation  data  at  individual  stressing  rates  were  negative.  The  n  parameter 
calculated  from  the  1000  gram  indentation  in  100%  RH  varied  from  4  at  the  120  MPa/sec 
stressing  rate  to  21  at  the  0.15  MPa/sec  stressing  rate.  While  the  latter  n  parameter  is  close  to 
that  determined  by  the  median  value  technique,  the  correlation  coefficient  was  0.56.  The  n 
parameter  values  for  the  50%  RH  data  varied  from  15  to  23.  These  are  considerably  less  than 
the  value  determined  by  the  median  value  technique. 


The  regression  analysis  of  d)mamic  fatigue  data  from  all  of  the  stressing  rates,  for  a  given 
indentation  load  and  environment,  using  the  homologous  ratio  gave  parameters  equivalent  to 
using  the  median  values  technique.  The  dynamic  fatigue  homologous  ratio  data  is  presented  in 
Figures  10-12. 


The  homologous  dynamic  fatigue  parameters  determined  for  individual  stressing  rates  bore  no 
resemblance  to  the  values  determined  using  the  median  values.  The  basic  assumption  that  the 
homologous  technique  hinges  on  is  that  the  strength  distribution  of  inert  specimens  and 
specimens  failing  at  various  stressing  rates  describe  the  same  flaw  population:  i.e.  the  weakest 
specimen  at  a  given  stressing  rate  will  fail  from  the  same  sized  initial  flaw  as  that  of  the  weakest 
inert  specimen.  The  use  of  the  indentation  technique  decreases  the  scatter  in  both  inert  and 
dynamic  fatigue  data.  The  scatter  in  the  static  fatigue  time  to  failure  is  substantial  even  with  the 
Vickers  indentations.  The  poor  agreement  obtained  using  the  homologous  ratio  analysis  of  the 
dynamic  fatigue  data  may  be  attributed  to  the  small  amount  of  scatter  in  both  the  inert  and 
dynamic  fatigue  data  sets. 
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Figure  9.  Static  fatigue  data  analyzed  using  the  homologous  ratio  technique.  The  line 
through  the  data  is  a  best  fit  using  all  results  from  the  four  applied  stresses. 


Figure  10.  Dynamic  fatigue  data  (4.9  N  indentation  in  water)  analyzed  using 
homologous  ratio  technique. 
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LN  [RATE/in^3(l/(s«:*Mpa"'2)] 


Figure  11.  E>ynamic  fatigue  data  (9,8  N  indentation  in  water)  analy2:ed  using 
homologous  ratio  technique. 


Figure  12.  Dynamic  fadgue  data  (9.8  N  indenUtion  in  50%  RH)  analyzed  using 
homologous  ratio  technique. 
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